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Abstract 

The dynamics of a degree of freedom associated to an axial vector in contact with 
a heat bath is described by means of a probabihty distribution function obey- 
ing a Fokker-Planck equation. The equation is derived by using mesoscopic non- 
equihbrium thermodynamics and permits a formulation of a dynamical theory for 
the axial degree of freedom (orientation, polarization) and its associated order pa- 
rameter. The theory is used to describe dielectric relaxation in homogeneous and 
non- homogeneous systems in the presence of strong electric fields. In the homoge- 
neous case, we obtain the dependence of the relaxation time on the external field as 
observed in experiments. In the non-homogeneous case, our model accounts for the 
two observed maxima of the dielectric loss giving a good quantitative description 
of experimental data at all frequencies, especially for systems with low molecular 
mass. 

Key words: Mesoscopic entropy, Fokker-Planck equations, axial degrees of 
freedom, dielectric relaxation. 



1 Introduction 



Orientational dynamics of a nematic director or of a dipole in an applied 
electric field plays a crucial role in numerous modern technologies. For exam- 
ple, the reorientation of the director vector in nematic liquid-crystals changes 
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the electro-optical properties of these materials, making of central interest 
the response of these systems to constant and time varying electrical fields 
[IJ. Liquid crystals, polymer solutions and dielectrics are systems showing 
strong anisotropics in their physical properties such as viscosities, transla- 
tional and rotational diffusion coefficients, and dielectric or magnetic suscep- 
tibilities p!][2][3]Hll5]l6|71l8]l9]IT0] . This fact together with the interaction of these 
systems with their environment, leads to a number of phenomena not ob- 
served in isotropic phases, like the dependence of the characteristic relaxation 
times of the reorientation of molecules on the intensity of an externally applied 
electrical field [TTIIT^ . 

From a fundamental viewpoint, the interaction with the surroundings and the 
orientational relaxation of the system can be modelled in terms of the dynam- 
ics of an axial degree of freedom n in contact with a heat bath [lll2] . Due to the 
importance of thermal fluctuations in the dynamics of the axial vector, in this 
article we begin our analysis by calculating the mesoscopic entropy production 
a of the system using the rules of mesoscopic nonequilibrium thermodynamics 
|13] . Then, from the entropy production and the probability conservation, we 
derive a Fokker-Planck equation for the probability density /(r, n, t) which 
depends on time t, position r, and the axial vector h. This Fokker-Planck 
equation allows us to formulate evolution equations for the average of the 
axial vector and for the associated order parameter in two important physi- 
cal situations: i) When the systems can be considered homogeneous and ii) 
when non-homogeneities are important. The results we have obtained can be 
compared with experiments [T^P^ and previous theories based on other ap- 

The paper is organized as follows. In Sec. 2 we calculate the mesoscopic entropy 
production, derive the Fokker-Planck equation for the probability density and 
obtain the evolution equations for the axial degree of freedom and its associ- 
ated order parameter under the assumption of the system's homogeneity. In 
Sec. 3 we apply our results to the case of dielectric relaxation and in Sec. 4 
we compare with experiments and with other theories. In sec. 5 we analyze 
the relaxation dynamics in non-homogeneous systems by taking into account 
the diffusion of the polarization. Finally, in the last section we summarize our 
main results. 



2 Mesoscopic entropy and Fokker-Planck dynamics 

Consider a system of particles in contact with a heat bath at constant tem- 
perature T and in the presence of an external potential U . The state of the 
system is characterized by the position vector r and the axial degree of freedom 
h which may be associated to an orientation or a polarization vector. 
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The mesoscopic dynamics of the system can be described in terms of the 
single-particle probability density f{f,n,t) satisfying the normalization con- 
dition / fdhdr = 1. This condition imposes that / must obey the probability 
conservation equation [T7] 



where 7? = ri x ^ is the rotational operator and V is the velocity conjugated 
to the position vector r. Eq. ([1]) can be obtained by taking into account the 
corresponding boundary and periodic conditions over f and n after applying 
the triple product identity and considering that the change in time of the axial 
vector n is perpendicular to fi: n = x n, with VL the angular velocity [2]. 

The explicit expressions for the unknown probability currents fV and fVt 
in (r*, n)-space, can be found by calculating the entropy production a of the 
system during its evolution in time. To this end, we will follow the rules of 
mesoscopic nonequilibrium thermodynamics which makes use of the Gibbs 
entropy postulate (GEP) p^llTBj 

ps{t) = -kB J f\n (^j-j dn + pSeq, (2) 



where ks is the Boltzmann constant, p = m J fdn is the mass density with m 
the mass of a particle and Seq is the equilibrium entropy, feq is the probability 
density characterizing the equilibrium state of the system, given by [T9|20|21] 

feq (r, n) = /oe-^/'^-^, (3) 



where, U = U {f, n) is a potential energy field per unit of mass and /o is a 
normalization constant. 

The GEP formulated in ([2]) is compatible with the Gibbs equation [18] 

T5s = 6e + p6 (p^^) - m J j26 (^^^ dndf, (4) 



where s, e, p and /i(r, h, t) are the non-equilibrium entropy, energy, pressure 
and chemical potential per unit mass of the system. Taking the variation of 
Eq. ([2]) and comparing with (jll) keeping constant e and p, we may identify the 
nonequilibrium chemical potential 

P = Peq + \n^ + U, (5) 

m to 
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where jleq is the chemical potential in the equilibrium state. The balance equa- 
tion for the non-equilibrium entropy follows from the Gibbs equation (with 
e and p constant) by taking the time derivative, integrating by parts and us- 
ing the corresponding boundary and periodic conditions for the currents. We 
obtain 



m 



/f|-(/^)-^'?A-(/«) 



dhdf. 



(6) 



Eq. (E]) contains two contributions of the type fiux-force pair, each one coming 
from the diffusion process in r- and n-space, respectively. Following nonequi- 
librium thermodynamics rules, linear phenomenological relationships between 
currents and thermodynamic forces can be established in the form p^p8] 



dp, 

or 

or 



(7) 



are the (tensorial) Onsager coefficients which are, 
in general, functions of r, n, T and p, [T31IT8] . However, for simplicity we will 
assume here that the Onsager tensors have constant components. Since n is an 
axial vector and r a polar vector, the Onsager- Casimir relationships impose 
the relation L„r = -L™, [ T8|l22] . 



Substituting Eqs. 
equation 



and (IE]) into ([T]), we finally arrive at the Fokker-Planck 



dt 



d_ 

df 



m 
m 



' df mf dU^ 
^ df ksT df ^ 



(9) 



Eq. ([9]) constitutes the basis of the mesoscopic description of the system and 
enables one to derive the evolution equation of the moments of the proba- 
bility density. A similar equation has been obtained by using hydrodynamic 
methods in Ref. [2] and from a Langevin approach in the case of a mean field 
approximation describing inhomogeneous nematics [23]. The present deriva- 
tion allows the generalization of Eq. ([9]) to cases in which inertia effects could 
be important, E 
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2.1 Relaxation equations for the average value of the axial degree of freedom 
and the order parameter 



In order to derive the relaxation equations for the average axial vector 
and the order parameter S, we will assume that the system is homogeneous, 
that is, the interaction potential of the particles with an external force only 
depends on n, f/ = U{h). In this case, the translational part of Eq. ([9]) can 
be integrated over ri, yielding a free diffusion equation for the center of mass 
of the particles from which the evolution equation for the second moment of 
position, (f^) = / f^f{r,n,t)dndr, is 

= 2Tr[Dt,]. (10) 



Here, Tr[Dtr] is the trace of the translational diffusion tensor I^tr assumed 
constant for simplicity and which is related to the Onsager coefficient L^^ 
through the expression = {ksT /myLrr-, |2]- Eq- f lTO|) implies that the center 
of mass of the particles performs Brownian motion with second moment (r^) = 
2Tr[Dt^]t. In the case of sohd systems at sufficiently low temperatures, the 
translational diffusion coefficient can be taken as zero. 

Thus, by integrating Eq. ([9]) over r we obtain the following Fokker-Planck 
equation for the reduced probability distribution g {h, t) = J fdr 



dg_ 
dt 



knT, 



m 



- mg ^ ~ 
Kg + TT^nU 



In the present case, the Onsager coefficient L„„ is related with the matrix 
of rotational relaxation times 7 through L„„ = which defines the 

rotational diffusion coefficient. From Eq. (fTT]) . we can obtain the relaxation 
equations for the first and second-order moments of g {n, t) corresponding to 
the average axial vector N and the order parameter S, defined by [2j 

N{t) = f ngdn and §(t) = (S)(t) = (1/2) / (3nn - l)gdh, (12) 



where 1 is the unit tensor and S = (l/2)(3?7,r;. — 1) is the microscopic order 
parameter. It is convenient to mention that when n represents the orientation 
vector of a molecule, then S is the orientation tensor [2^1] that is related 
to the well known Onsager order parameter characterizing the order- disorder 
transition in suspensions of rod-like molecules [3]. For polar systems, n and § 
are related to the polar and quadrupolar terms of the corresponding multipolar 
expansion. 
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Now, taking the time derivative of equation ( !T2|) and using ( !TT|) . after inte- 
grating by parts assuming periodic conditions for the probabihty current, one 
obtains the evolution equation 



ON 
~dt 



3kBT- \ ^ ^' dh/^ 3kBT \ dn ^ 



(13) 



where M(S) is defined by 

M(S) = [7-1 ■ S + (j-' : S) 1 - • S 



(14) 



In Eq. f|T3l) . the matrix of relaxation times is given by 7"^ = Tr(7~^)l — 
7~^, implying that, in the general case, the characteristic relaxation times are 
combinations of the elements of the rotational diffusion tensor. The first term 
on the right-hand side of Eq. ( fT3l) comes from the diffusion term of Eq. ( fTTl) , 
whereas the other terms are due to the external force. Eq. ( [T3l) is expected to 
be valid for characteristic values of the interaction energy, U, of the order of 
the thermal energy /c^T. 



Following a similar procedure, Eqs. ( ITTll and ( !T2l) yield the evolution equation 
for § 



as 



-2M-(S)-2(r'.Sr + ^(M(Q).f~" 





on 



(15) 



where the superscript s stands for the symmetric part of a tensor and M(S) 
has the same form as in Eq. (IT^ but in (ITSl) depends on the averaged order 
parameter S. In Eq. (fT5|) we have introduced the average tensor Q defined by 



Q(t) = (Q)(t) = (1/2) / [5nnh - {nlY]gdn 



(16) 



and the tensor 



M(Q) = 7 



-1-Q+(7-^:Q)1-7-'-Q 



(17) 



Here, Q is the microscopic octupolar tensor and (^l)^^;. = ni6jk + nj6ki + nkSij. 
In Eq. (ITSll . the first three terms come from the diffusion term in Eq. ( ITTl) 



6 



whereas the others are due to interactions with the external field. Eqs. (fT3!l 
and (fTSll constitute the first two evolution equations of the complete hierarchy 
of moments of the distribution function g. 



3 Dipole and quadrupole relaxation in a homogeneous system 



The relaxation of the polarization vector in dielectrics can be described by 
means of equations (fT3|) and (fT5l) since h can be interpreted as the dipole 
moment of a molecule through the expression: n = n~^p where n is the mag- 
nitude of the dipole moment of the molecule. In this form, the polarization P 
and the cuadrupolar moment Sp per unit volume are defined by 



P = Nj^N and §p = A^>^ (S) 



(18) 



where the tensor Sp is the quadrupole moment induced by thermal fiuctua- 
tions. In first approximation, the potential energy per unit mass of a polar 
molecule interacting with an external electric field E is given by [2j 



p- E 



m 



(19) 



From Eqs. ( IT^ . (113]) and (ITS]) , it follows that the evolution equation for P is 



dP 
'dt 



P + J 7 ■ M(§.) ■ E - ^^E 

3NkBT- ^ ^' SkeT 



(20) 



where M(Sp) is defined according to Eq. (JHj). The evolution equation for Sp 
takes the form 



'dt 



-2M'(§p) -2 (7-^ -Sp)' 
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5NkBT 



■E 



3Nj? 
'bkeT 

'bksT 



{2 [(r'-^)^]' + [r' {p-^)Y} 

{[{i'-e)-p]i + [{i'-p) e]'} 



(21) 



where Qp = N^n^Q with M(Qp) defined according to Eq. i^. The con- 
tribution of higher order moments will be neglected for convenience. These 
equations generalize the Debye theory of dielectric relaxation by taking into 
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account the coupling between the dipolar and quadrupolar moments natu- 
rally contained in the distribution function g. Debye's theory is recovered if 
quadrupolar contributions are neglected in Eqs. ([20]) and (PTi) . 

A closed set of equations for P and E>p can be obtained by approximating 
Qp by its equilibrium average, which can be calculated using the equilibrium 
distribution function 

ge, = Z-^e^'^l^-^. (22) 



Eq. ([22]) is the stationary solution of ^ and Z"^ = J e"^^ '^'^'^ dp is the 
partition function of the system. Notice that, when the externally applied field 
is weak, it is usual to take the so-called homogeneous isotropic approximation 
of ([22]), given by [2] 

^ (23) 



We will use Eqs. (1221) and ( l23l) in order to complete the set of equations (1201) 
and (l2T]) and to make comparisons with previous theories. 



3.1 Dielectric relaxation in the homogeneous-isotropic approximation 



By using Eqs. f[20]) and f[2T]) . and the isotropic distribution function f[23]) . the 
dielectric response of a polar system will be analyzed in terms of the complex 
dielectric susceptibihty x- 

To proceed, we will now assume that the external electric field is time depen- 
dent and given by the sum of a strong constant electric field Eq applied along 
the z axis, and a weak AC probe field Ei{t) = i^ie"^"* applied in an arbitrary 
direction, that is: E = Eq + Ei{t). In addition, we will assume for simplicity 
that dE>p/ (9t ~ and solve Eqs. (!20!) and (pT!) for P. This quasi-stationary as- 
sumption leads to a considerable simplification of the obtained formulas and 
do not modifies considerably the results of our model. Thus, since the average 
of the octupolar tensor ( TTBi) with ([25]) vanishes ([Qpjijfc = 0), from Eq. we 
obtain the following relation for Sp 



'(Sp) + 2 7 



{2 [(r^ ■ m) p\ 



p-m] 



■ m) ■ p] 1} . (24) 



8 



From Eqs. (12^ and (120|) . it is possible to obtain two equations for the parallel 
P\\{t) and perpendicular P±{t) components of the polarization vector defined 
with respect to the constant field Eq. Hence, using the condition Ei « Eq 
in the second term of Eq. (1201) and keeping dominant terms, the equation for 
P\\{t) becomes 



gP||(t) 
dt 



-Til 



(25) 



where E\\ = Eq + Ef with Ej the component of Ei paralel to Eq, and we have 
assumed that the matrix of relaxation times 7 is diagonal with components 7^, 
7_L and 7||. By definition, 7 and 7 have the same symmetries. The relaxation 
time T\\ and the zero frequency susceptibility x\\ iii Eq. ( 125|) are given by 



Til = 75 and Y|| = To . 26 

i + i(4-r) ''I' i + i(4-r) 



Here, we have introduced F = 7_l/7|| to simplify the notation and defined 
Xo = N^^/'itokBT with eo = 8.854 x 10~^^F/m, [25]. It is important to notice 
the fact that the constant electric field introduces corrections to r|| and x\\ 
through the dimensionless parameter ^ = fiEQ/ksT. 

Using Fourier transform techniques, from Eq. (12 5 p one obtains 

P|,(t) = eoX||^o + eoX||K)4'W, (27) 



where X\\{^o) = + iu}oT\\)~^- The evolution equation for the perpendicular 
component of the polarization vector can be calculated in similar form. One 
obtains 



dt 



-7± 



P^it) - eoXoE[-^\t) 



(28) 



where e['^'^ is the component of Ei perpendicular to Eq. Following the same 
procedure to obtain f l26l) . we find that: Pj_ 
case, r_L and x± take the form 



eoXL{l + ii^or±)-^E[^\t). In this 



7±/(i + r) 

10 (i+r) 



and 



X± 



1 + 



«2 (2-r) 
10 (i+r) 



(29) 



These results yield the following expressions for the real and imaginary parts 
of the parallel and perpendicular components of the complex dielectric sus- 
ceptibility tensor 
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xiM = ^^ + xr and xj[M = ^^^, (30) 
xiK) = -^ + xr and ^'A^o) = ^f^. (31) 

where xf and Xj? are the high frequency susceptibihties of the system. From 
Eqs. fl26l) and fl29l) . it follows that the quadrupole moment induced by thermal 
flucutations introduces field dependent corrections in the susceptibility and the 
relaxation times, as observed in experiments [T2] . 



We have also calculated the expressions for the parallel and perpendicular 
susceptibilities of the system in the non-stationary case. The expressions for 
the components of the polarization vector can only be given in Fourier space 
and, in the isotropic approximation they are given by 

^,TLS ns 

Pr(w) = Co Eu(uj) and Prico) = Co — EAuj), (32) 



where the superscript ns indicates the non-stationary case. The relaxation 
times Tu"'' and T]f are now given by 



^ns _ 7±/2 7±/(l + r) 

'II ~ , , ^2 (4_r) anu /_|_ — ^ c2 (2-r) ' v'-''-'/ 

30 [l+iuiT^/efl + To (l+r)[l+(Lrx/6)2J 



whereas the susceptibilities x^ and x"'' in this non- stationary case are defined 
by 

^ + 30 [l+(wrx/6)"J ^ + 10 {l+r)[l+{ujT^/6f] 

Finally, the correction factors c\\ and c± due to the time derivative in Eq. fl2Tl) 
are respectively given by 



(4-r 



90 [l + {uT^/Qf 

e (2-r) 



60 l + {ujTj6y 



(35) 



In the following section, the comparison with experiments using the complete 
set of non-stationary relations for the susceptibilities, Eqs. ( 132|) - (!35|) . or their 
non-homogeneous versions (see Sec. 5), will show that the quantitative descrip- 
tion is very similar to the one given by means of the stationary expressions, 
that are simpler. 
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3.2 Relaxation in an homogeneous non-isotropic approximation 



Eqs. fl26|) and f l29|) for the relaxation times and the susceptibihties have the 
correct quahtative dependence on the external field Eq as observed in exper- 
iments, but they fail quantitatively. This is due to the fact that the isotropic 
approximation is not adequate when strong electrical fields are applied to the 
system. 

Therefore, the quantitative description can be improved by taking the average 
f ll6p in Eqs. f l20p and f l2ip with the non-isotropic probability distribution fl22p . 
In this case, the components of the tensor Qp do not vanish and more compli- 
cated formulas are obtained for Sp and P. Thus, following a procedure similar 
to that of subsection 3.1, it is possible to derive the following expressions for 
the zero frequency parallel and perpendicular susceptibilities 



XII 



xo{i + f(3-2r) [{i + 'i)m-i]} 



i + lo(4-r) 



Xo 






{m - 2r ' 




m - i 


}] 




.-^ 10 (1+r). 





(36) 



Here, the tilde stands for the non-isotropic case. The expression for the re- 
laxation times do not change and are given by the corresponding expressions 
in (j26l) and (j29ll . If we approximate our description by decoupling Eqs. (!20l) 



and (I2T1) by taking the average of Sp with ( l22l) . then the expressions for the 
components of the susceptibility reduce to 

^xom] 



~bh 

X\\ 



and 



x't 



Xo 




, (37) 



where the parallel component xjl'* is equal to the one reported in Ref. [12]. In 
this approximation, the parallel and perpendicular relaxation times become: 
ry =7^/2 andr^ = 7^/[l + r]. 

In this section we have shown that the field dependent corrections of the 
relaxation times and susceptibilities of a dielectric are due to the coupling 
between the evolution equations for the dipolar, quadrupolar and, in general, 
higher order moments of the distribution. The results obtained agree with 
previous theories (see, for instance, [16]) and, when ~ 0, reduce to the ones 
derived by means of Debye's theory. 

When considering the non-stationary dynamics of the quadrupole, we obtain 
that in the present case the susceptibilities entering in the corresponding Eqs. 
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Fig. 1. Comparison of the experimental data of the critical frequency fc{E) (sym- 
bols) taken from Ref. [12] with the theoretical results for a strong DC field 
(lines) for a dilute solution of PBLG molecules with molecular weights of 1.1 
xlO^(romboids), 2.6 xlO^(circles) and 4.6 x 10^ (squares). 

(l32l) are now given by 



{1 + 



xT = Xo 



{1 + 

xr=xo^ — 



g (3-2r) 

5 (1— iwTx/C) 



(i + P)iKl-f]} 



1 + (4-r) 

30 \l+(u^r^/&f\ 



(2-r) 



5 (l+r)(l-iwTx/6) 



{^[e]-r[(i + f|)L[e]-f]}} 



1 + 



$2 [(2-r)/(i+r)] 

10 [l+(a;Tx/6)2] 



(38) 



The relaxation times and the correction factors cy, c_l are still given by Eqs. 
fl33l) and fl35|) . In the case when there is no coupling we have cy = c_l = 1, and 
therefore the relaxation times and the susceptibilities are given reduce to the 
form Eq. fl37|) . as expected. 



4 Comparison with experiments 



In this section, we use Eq. fl26|) for the relaxation time to estimate the value of 
the dipolar moment of PBLG molecules with different molecular masses [12j. 
In similar form, the susceptibilities given in f l36p are used to fit the experi- 
mental data for the dielectric loss of czs-PI polymers in solution reported in 
Ref. [26]. 



4-1 Determination of the dipolar moment using the field dependent relaxation 
time 



The magnitude of the dipolar moment /i of a molecule can be estimated by 
using the critical frequency /c, defined as the value of the frequency at which 
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M„(xl05) 


Debye 


BH-1 


BH-2 


WKC 


UUman 


FP 


4.6 


8.7 


7.1 


26 


23.3 


22.3 


23.2 


2.6 


5.6 


4.6 




14.7 




13.9 


1.1 


2.8 


2.3 




4.8 




4.0 



Table 1 

Values of the dipolar moment fi (lO^^^C m) of a dilute solution of PBLG polar 
molecules having different molecular masses: = 1.1, 2.6 and 4.6 x 10^ The first 
line denotes different models used for comparison (see text for details). The last 
column contains the results obtained with our model. 

the maximum of the dielectric loss occurs. Since the value of fc depends on 
the magnitude of the applied field, see Fig. 1, then fi can be determined by 
using the relation [8.26j 

/, = (27rr||)-\ (39) 

The open symbols in Fig. 1 show isothermal (298 K) data of the critical 
frequency fc, reported from experiments with dilute solutions of PBLG having 
different molecular masses M^,, [I2]. The data were obtained by applying a 
strong DC electric field (0 — 9 x 10^ V/m) on the system and can be fitted 
using Eq. (126|) by assuming that the energy of a dipole in the DC field is of 
the order of the thermal energy (^ ^ 1), and that 7_l ~ 7|| = 7. In this case 
the relaxation time takes the form ry = (7/2) /[I + (,^^/10)]. 

Table 1 shows a comparison of the values of fi obtained by using different 
theories. The values estimated with the Debye and the Block and Hayes model 
1 (BH-1), were obtained in the weak- field approximation. For the system with 
the largest molecular mass in the strong-field approximation, the magnitude 
of the dipolar moment obtained with the Block and Hayes model 2 (BH-2) is 
very similar to those obtained by means of a theory based on the Langevin 
equation (WKC) |16], and a classical result reported by Ullman in Ref. [15]. 
The last column in Table 1 shows that the results obtained with our model 
(FP) in the non-isotropic (strong field) approximation are in good agreement 
with the experimental results. 

4-2 Normal and a relaxations in type-C polymers 

Highly concentrated solutions of czs-polyisoprene (cis-Pl) show a dielectric 
relaxation that involves both normal and a-relaxations characterized by two 
absorption peaks of the imaginary part of the dielectric susceptibility; see Fig. 
[2] and Ref. [26] . The first maximum of the dielectric loss is usually attributed 
to the dipolar moment aligned with the chain contour, and is associated to 
the normal relaxation characteristic of Type-A polymers. The second maxi- 
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-1 0123456789 10 11 -1 0123456789 10 11 



log^o(m) log^^(») 



0.04 




-1 0123456789 10 11 



log^„(«) 

Fig. 2. Experimental data (symbols) [26j of the dielectric loss of highly concen- 
trated solution of cis-PI chains at 273 K as a function of frequency for different 
mass weights The average molecular mass in kg-mol~^ are: a) 2; b) 14 and c) 
32. The dotted lines correspond to the best fit obtained with Eqs. (pOj) and (j3T|) 
under the assumption of an homogeneous system. The solid lines correspond to the 
best fit using Eqs. ()43p and ()44p for the non-homogeneous approximation where 
we have used dispersion relations of the form k'^{uj) = (wq/w*)" for the parallel 
and perpendicular components of the translational diffusion coefficient. Here, a ^ 1 
and uj* is a characteristic frequency associated with the viscoelastic nature of the 
medium (see details in Sec. 5). The dashed lines correspond to the non-stationary 
and non- homogeneous approximation given by Eqs. (I45p and (j46p . 



mum is related to the dipolar moment perpendicular to the chain contour and 
characterizes Type-B polymers. The cis-Pl presents both types of relaxations 
and therefore is classified as a Type-C polymer. 

For this system, the polarization p can be interpreted as a local degree of 
freedom in contact with a heat bath, and therefore its dielectric relaxation 
can be described by using Eqs. (!30|) - (!36|) only after assuming that the system 
is homogeneous (the translational diffusion of the polarization is neglected). 

The comparison between experiments (symbols) and theory (lines) shown in 
Fig. Ill was performed by taking the total imaginary part of the susceptibility 
as the sum of the parallel X\\{^o) and perpendicular x±{^o) contributions given 
in ( l30l) and ( 13T|) . respectively. The low frequency maximum of x" was used 
to obtain the values of ^, Xo and ry as was done in section 4.1. Then, after 
determining the characteristic frequency (tJ^) at which the high frequency 
maximum appears, Eqs. (130|) and ( 13T|) were used to fit the whole set of experi- 
mental data. From the dotted line in Fig. it is clear that our model gives an 
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Polymer 


XO 




7|| [s] 


7± [s] 


^11 [^1 


r±[s] 


PI02 


0.13 


7.48 


1.66x10-5 


8.25x10"* 


1.96x10-*^ 


6.77x10-^ 


PI14 


0.13 


7.98 


6.18x10"^ 


1.69x10-^ 


6.51x10"* 


1.23x10-* 


PI32 


0.12 


7.07 


1.24x10-1 


1.69x10-7 


1.62x10-2 


1.54x10-* 



Table 2 

Values of the parameters obtained with Eqs. ()30p and ()3ip in the homogeneous 
approximation (see dotted lines in Fig. [2]). xo and ^ are dimensionless quantities. 

amplitude of the second maximum which coincides well with the experimental 
one for the three molecular masses of the cis-Pl considered. Table 2 contains 
the values of the parameters used to fit the data. Notice that, in this case, 
the interaction energy per molecule between the system and the external con- 
stant field is about five times larger than the thermal energy. Since only two 
relaxation modes were considered, the theoretical curves are narrower than 
the experimental ones, as expected. In the following section we will show that, 
to obtain broader curves, one has to consider the translational diffusion of 
polarization. 



5 Diffusion effects in dielectric relaxation 



The previous description of dielectric relaxation in cis-Fl solutions can be 
improved by considering the effects of non-homogeneities of the dipolar and the 
quadrupolar moments. Therefore, it is relevant to consider the local definitions 



P(r, t) = N I pf{p, r, t)dp; §p(r, t) = — 



i3pp-fxH)f{p,r,t)dp, (40) 



where we have assumed that the distribution function depends on the polar- 
ization vector 

Similarly as in section 2, the evolution equations for P{f,t) and E>p{f,t) can 
be calculated by taking the time derivative of Eqs. (HU]1 . substituting in 
the resulting expression and integrating by parts assuming the corresponding 
boundary conditions. The resulting diffusion equation for P{f, t) is 



dP 
'dt 



Hr 



: VV) P - 7 



p + 



1 



whereas for §p(r, t) the diffusion equation is 



(41) 



d8p 
'dt 



Hr 



: VV) §p 



6 

bNkeT 



M(Qp) -E]' - 2M*(§p) - 2 (7-^ • 
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Polymer 


XO 




7|| (s) 


7±(s) 


a 


P 


a;||(Hz) 


wx{Hz) 


PI02 


1.48 


2.51 


4.10x10-*^ 


5.75x10-* 


1.03 


0.65 


2.03x10-^ 


3.39x10-^ 


PI14 


8.14 


2.15 


2.93x10-* 


1.00x10-* 


0.79 


0.71 


2.19x10-^ 


3.38x10-^* 


PI32 


3.69 


1.87 


9.79x10-2 


3.73x10-* 


0.80 


0.65 


1.10x10-2 


1.12x10-^* 



Table 3 

The numerical values are obtained with nonlinear regression with Eqs. (j43p and 
in the non-homogeneous approximation (dashed lines in Fig. [2]). xO) ol and (3 are 
dimensionless quantities. The values of the relaxation times r can be obtained by 
using Eqs. (l26]l and (f29l) . 



3iV 



{2 [{r' ■ e) P 

7 



+ 



7-1 [P-E 



--'■E] -P 



7 



(42) 



The first term on the right-hand side of equations ( 14T1) and ( l42l) accounts for 
the translational diffusion, V = d/df is the gradient operator and we have 
assumed that the interaction potential only depends on the externally applied 
electric field: U = —p- Eq. Eqs. ( HTj) and P2|) are coupled through the term 
containing the tensor M(§p). 

The expressions for the parallel and perpendicular contributions of the dielec- 
tric susceptibility tensor can be obtained by following a procedure similar to 
that of Sec. 3. That is, by substituting the components of §>p obtained from 
Eq. (1421) into (|4T1) . in the quasi- stationary approximation one finds 



X\\{k,uJo) 
X±{k, Wo) 



X\\ 



x± 

1 + Dtr±k^r± + iUor_L 



(43) 
(44) 



The expressions for the relaxation times ry and are given through Eqs. 
and f l29|) . The relation for the parallel and perpendicular components of the 
susceptibility x\\ X± are still given by f p6|) . 

Comparison with experiments can be performed by using the previous expres- 
sions and assuming that the parallel and perpendicular components of the 
diffusion tensor are of the form: Dfruk"^ = (cjo/i^ii)" and Dtri.k'^ = {uJo/oj_i_Y . 
This dependence on frequency has been suggested in the hterature and is re- 
lated to the viscoelastic nature of the system [27] , and discussed in relation to 
the problem of liquid-crystal like relaxing amphiphilic systems such as lipid 
membrane. In these last systems, it was shown that the non-Debye relaxation 
manifests a memory feeling submesoscopic dynamics that may arise due to 
a reduced system's dimensionality that can be accessed phenomenologically 
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In the non- stationary case, the expressions obtained for the susceptibihties are 



~ns 
~ns 

where all the field dependent factors (r"'^, T]f, cy, c±, X\\^ and xT) are given 
by Eqs. ([33]), §E) and (|38|). 

The solid line in Fig. |2] shows the theoretical results in the non-homogeneous 
case. The effect of diffusion is to make the peaks broader due to the power 
law relation between the wave number and the frquency. They also introduce 
two relaxation times associated to u^^ and wj^. Fig. 2 shows that our ap- 
proach is more appropriate for describing systems made up of particles with 
low molecular mass, although both maxima are well described even for sys- 
tems of particles having large molecular masses. Table 3 shows the values of 
the parameters obtained in the present approximation after using a nonlin- 
ear regression of Eqs. (H3|) and (144|) to the experimental data. Calculations 
in progress indicate that the description offered here could be improved by 
incorporating the effects of mass dispersion. 

The results obtained in this section suggest that the coupled relaxation equa- 
tions for the dipolar P and quadrupolar Sp fields constitute a minimal model 
to describe the dielectric relaxation of solutions of Type-C polymers. 



6 Conclusions 



We have formulated a mesoscopic model describing the dynamics of systems 
characterized by an axial degree of freedom in contact with a heat bath. The 
description of the dynamics of this degree of freedom is performed by means 
of a Fokker-Planck equation for the nonequilibrium probability distribution 
function. The equation, derived by calculating the mesoscopic entropy pro- 
duction, allows one to formulate a nonequilibrium description of the system 
in terms of relaxation equations for the multipoles, related to the moments of 
the distribution function. 

Applications of the theory to describe dielectric relaxation in homogeneous and 
non-homogeneous systems in different approximations allowed us to compare 
with previous theories and better elucidate the new contributions of our model. 
In the case of dilute polymer suspensions, we found that our theory accounts 
for the dependence of the relaxation time on the strong applied electrical field 
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and gives values for the dipolar moment of the particles which agree with 
previous ones. 

We used our model to describe the normal and a relaxations in Type-C poly- 
mer solutions. Our results clearly show that, in order to correctly account for 
the amplitude of the second maximum of the dielectric loss (associated to the 
a relaxation), the minimal model has to take into account the coupled diffu- 
sion equations for the dipolar and quadrupolar fields. The anisotropic diffusion 
of these quantities introduces new characteristic times leading to a better de- 
scription of the experiments. In particular, we found that our theory is very 
useful in describing dielectric relaxation of solutions of molecules having low 
molecular masses. 

The general model describing the dynamics of systems characterized by axial 
degrees of freedom presented here can be generalized and applied to other 
systems such as suspensions of rod-like particles or liquid crystals through an 
adequate choice of the interaction energy. This will be done in future work. 
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